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Abstract. We present a general theory of fractal transformations and show how it leads to new 
type of method for filtering and transforming digital images. This work substantially generalizes 
earlier work on fractal tops. The approach involves fractal geometry, chaotic dynamics, and an in- 
terplay between discrete and continuous representations. The underlying mathematics is established 
and applications to digital imaging are described and exemplified. 
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1. Introduction. Fractal transformations are mappings between pairs of at- 
tractors of iterated function systems. They are denned with the aid of code space 
structures, and can be quite simple to handle and compute. They can be applied to 
digital images when the attractors are rectangular subsets of K 2 . They are termed 
"fractal" because they can change the box-counting, Hausdorff, and other dimensions 
of sets and measures upon which they act. In this paper we substantially generalize 
and develop the theory and we illustrate how it may be applied to digital imaging. 
Previous work was restricted to fractal transformations defined using fractal tops. 

Fractal tops were introduced in [2] and further developed in [5J H3 [7J QT] . The 
main idea is this: given an iterated function system with a coding map and an at- 
tractor, a section of the coding map, called a tops function, can be defined using 
the "top" addresses of points on the attractor. Given two iterated function systems 
each with an attractor, a coding map, and a common code space, a mapping from 
one attractor to the other can be constructed by composing the tops function, for 
the first iterated function system, with the coding map for the second system. Under 
various conditions the composed map, from one attractor to the other, is continuous 
or a homeomorphism. In the cases of affine and projective iterated function systems, 
practical methods based on the chaos game algorithm [H] are feasible for the approx- 
imate digital computation of such transformations. Fractal tops have applications to 
information theory and to computer graphics. They have been applied to the produc- 
tion of artwork, as discussed for example in [4] , and to real-time image synthesis |18) . 
In the present paper we extend the theory and applications. 

Much of the material in this paper is new. The underlying new idea is that diverse 
sections of a coding map may be defined quite generally, but specifically enough to be 
useful, by associating certain dynamical systems with the iterated function system. 
These sections provide novel collections of fractal transformations; by their means 
we generalize the theory and applications of fractal tops. We establish properties of 
fractal transformations, including conditions under which they are continuous. The 
properties are illustrated by examples related to digital imaging. 

A notable result, Theorem |5.3[ states the existence of nontrivial fractal homeo- 
morphisms between attractors of some affine overlapping iterated function systems. 
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Fig. 1.1. Lena after application of a fractal homeomorphism. See Example \5.5\ 



The proof explains how to construct them. An example of one of these new homeo- 
morphisms, applied to a picture of Lena, is illustrated in Figure [L~T| 

In Section [2] we review briefly the key definitions and results concerning point- 
fibred iterated function systems on compact Hausdorff spaces. Since this material is 
not well-known, it is of independent interest. The main result is Theorem |2.4| This 
can be viewed as a restatement of some ideas in [13] : it describes the relationship 
between the coding map and the attractor of a point-fibred iterated function system. 

In Section [3] we define, and establish some general properties of, fractal trans- 
formations constructed using sections of coding maps. In Theorem |3.2| we present 
some general properties of coding maps. Then we use coding maps to define fractal 
transformations and, in Theorem 3.4 we provide sufficient conditions for a fractal 
transformation to be continuous or homcomorphic. 

In Section [4] we define two different types of section of a coding map: (i) with the 



aid of a masked dynamical system; and (ii) with the aid of fractal tops. Theorem 4.3 



establishes the connection between the masked dynamical system and a masked section 
of the coding map. Theorem 4.5 includes a statement concerning the relationship 
between the masked section to the coding map and the shift map. Here we also 
establish the relationship between fractal tops and masked systems. A key result, 
Theorem |4.7| gives a condition under which the ranges of different masked sections 
intersect in a set of measure zero. This enables the approximate storage of multiple 



images in a single image, as illustrated in Figure 5.7 



In Section[5]we apply and illustrate the theoretical structures of Sections 2,3, and 
4, in the context of digital imaging. Our goal is to illustate the diversity of imaging 
techniques that are made feasible by our techniques, and to suggest that fractal trans- 



formations have a potentially valuable role to play in digital imaging. In Section 5.1 
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we illustrate how fractal transformations may be applied to image synthesis, that is to 
making artificial interesting and even beautiful pictures. Specifically we explain how 



the technique of color-stealing [5] extends to masked systems. In Section 5.2 we apply 



fractal homeomorphisms, using Theorem 3.4 to transform digital images, for image 



beautification, roughening, and special effects; in particular, we present and illustrate 
Theorem |5.3| which extends the set of known affine fractal homeomophisms. In Sec- 
tion |5.3| we consider the idea of composing a fractal transformation, discretization, 
and the inverse of the transformation to make idempotent image filters. In Section[5.4| 



we apply Theorem 4.7 to the approximate storage or encryption of multiple images 



in a single image. In Section 5.5 we provide a second technique for combining several 
images in one: it combines invariant measures of a single iterated function system with 
several sets of probabilities, to make a single "encoded" image: approximations to the 
original images are revealed by the application of several fractal homeomorphisms. 

2. Point-fibred iterated function systems. Let X be a nonempty compact 
Hausdorff space, and let K (X) be the set of nonempty compact subsets of X. It is 
known that K(X) endowed with the Vietoris topology is a compact Hausdorff space, 
see for example [TH Theorem 2.3.5, p. 17]. This encompasses the well-known fact that 
if X is a compact metric space then K(X) endowed with the Hausdorff metric is a 
compact metric space. Let I = {1,2, ...N} be a finite index set with the discrete 
topology. Let {/j : X — > X\i 6 7} be a sequence of continuous functions. Following 

m, 

T:= (X-h,...f N ) 

is called an iterated function system over X. 

Following [13l Definition 4.1.4, p. 84] we define a map 

oo 

n : 7°° -> K{X), a ^ f| f ai o f„ 2 o • • • o f„ k (X) (2.1) 
fe=i 

for all sequences a — <j\ o^o?,-.- belonging to The map is well-defined because 
n(cr) is the intersection of a nested sequence of nonempty compact sets. The following 
definition is based on [13l Definition 4.3.6, p. 97]. 

Definition 2.1. Let J- :— (X; fx, ■ ■ ■ /at) be an iterated junction system over a 
compact Hausdorff space X . If II(er) is a singleton for all a € I°° then J- is said to 
be point-fibred, and the coding map of T is defined by 

tt : 7°° -> A, {vr(a)} = 11(a), 

where Ad X denotes the range ofir. 

Theorem |2.2[ due to Kieninger, plays a central role in this paper. It generalizes 
a classical result of Hutchinson [TU] that applies when X is a compact metric space 
and each / € T is a contraction. 

Theorem 2.2. Let I°° have the product topology. If J- is a point-fibred iterated 
function system on a compact Hausdorff space X then the coding map tt : I°° — > A is 
continuous. 

Proof. This follows from [T31 Proposition 4.3.22, p. 105]. □ 
We define 

T : K{X) -> K(X), B >-» \J f{B). 



4 



M. F. BARNSLEY, B. HARDING, K. IGUDESMAN 



By slight abuse of notation we use the same symbol T for the iterated function system, 
the maps that it comprises, and the latter function. We define F a = ix, the identity 
map on X, and F k = J- o J rk ~ 1 for k = 1, 2, ... . The following definition is a natural 
generalization of the notion of an attractor of a contractive iterated function system, 
see for example [15l definition on p. 1193 and Theorem 11.1, p. 1206], [13l p. 107], and 
also [8]. 

Definition 2.3. Let T be an iterated junction system on a compact Hausdorff 
space X. An attractor of J- is A G K(X) with these properties: (i) J^{A) = A; (ii) 
there exists an open set U C X such that A C U and 

lim F k {B) = A (2.2) 

k— f oo 

for all B C U with B G K(X). (The limit is with respect to the Vietoris topology on 



K(X).) The largest open set 6cl such that equation 2.2 holds for all B C B with 
B £ K(X) is called the basin of A. 

The relationship between coding maps and attractors is provided by noting that 

tt(ct) = lim f a \ k {a) where f a \ k := f ai o f a2 o ■ ■ ■ o / 

k— ^oo 

for a G X. This leads to the following theorem. 

Theorem 2.4. If J- is a point-fibred iterated function system on a compact 
Hausdorff space X then 

(i) T : K(X) — > K{X) has a unique fixed-point A G K(X), i.e. J 7 (A) = A; 

(ii) A is the unique attractor of T: 

(Hi) A is equal to the range of the coding map ir, namely 

A = 7r(7°°); 

(iv) the basin of A is X; 

(v) ifBe K{X) then {ir{a)} = linu^^ f a \ k (B) for all a G I°° . 

Proof. This follows from O Proposition 4.4.2, p. 107, see also Proposition 3.4.4, 
p.77]. □ 

The following remark tells us that if an iterated function system possesses an at- 
tractor then it is point-fibred when restricted to a certain neighborhood of the attrac- 
tor. Necessary and sufficient conditions for an iterated function system of projective 
transformations to possess an attractor are given in [S]. 

Remark 2.5. Let A be an attractor of an iterated function system T on a compact 
Hausdorff space X , and let B be the basin of A. Then, following [15], T defines an 
iterated closed relation r :— {(x,f(x)) : x G X, f G J-} C X x X, and A is an 
attractor of r. By \ 151 Theorem 7.2, p. 1193] there exists a compact neighborhood V 
of A such that A C Intx(V) and J- (V) C Intx(V). (Intx(S) denotes the interior 
of the set S G X in the subspace topology induced on S by X.) It follows that 
J-\v '■= (V; /i|vj f2\v, ■■•/iv|v) is a point-fibred iterated function system on a compact 
Hausdorff space. The attractor of J-\y is A. 

Definition 2.6. Let T be a point-fibred iterated function system on a compact 
Hausdorff space. The set I°° is called the code space of J 7 . A point a G I°° is called 
an address of it (a) G A. 

In the rest of this paper the underlying space X is a compact Hausdorff space. 
Also in the rest of this paper the symbols J r ,G,'H denote point-fibred iterated function 
systems on compact Hausdorff spaces. We will say that an iterated function system 
is injective when all of the maps that it comprises are injective. We will say that an 
iterated function system is open when all of the maps that it comprises are open. 
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3. Fractal transformations. Here we present a generalized theory of fractal 
transformations and establish some continuity properties. Fractal transformations are 
defined using sections of coding maps. We are concerned with continuity properties; 



for example, Theorem 3.4 (i) provides a sufficient condition for a fractal transformation 
to be continuous. 

Definition 3.1. Let ir : I°° — >• A be the coding map of T . A subset SI C I°° is 
called an address space for J- if Tt(Sl) = A and 7r|n : SI — > A C X is one-to-one. 
The corresponding map 

t : A Sl,x m- (Tr\ n y 1 (x), 

is called a section ofir. 

Theorem |3.2| summarises the properties of sections of ir. 

Theorem 3.2. Let T — (X; /i, f%, /jy) be a point-fibred iterated function 
system on a compact Hausdorff space, with attractor A, code space I°° , and coding 
map 7r : I°° — > A. If r : A — > SI is a section of tt then 

(i) t : A SI is bijective; 

(ii) t -1 : SI — » A is continuous; 

(Hi) iot = la, the identity map on A, and r o (7r|n) = i^n, the identity map on 



Si; 



(iv) if T is infective and fi(A)C]fj(A) = for all i,j £ I with i ^ j , then SI = 1° 

(v) if SI is closed then r : A — > SI is a homeomorphism; 

(vi) if A is connected and A is not a singleton, then r : A — > SI is not continuous. 



Proof, (i) By Definition 3.1 7r| o : SI — > A is bijective, so r = (7r|o) _1 : A — >• SI is 
bijective. 

(ii) By Theorem 2.2 ir : I°° -> A is continuous. It follows that r _1 = tt\q : SI — s- A 
is continuous. 

(iii) If x £ A then ir o t(x) — ir o (7t|q) 1 (x) C tt o 7r _1 (a;) = x = iA{%)- Also 
r o 7r|o = r o r _1 = lq. 

(iv) Suppose 17 ^ I°° . Then there are a, uj € 7°°, cr 7^ w, such that 7r(er) = 
7r(w). We show that this is impossible. If ir{a) = n(u)) then Theorem |2.4| (v) implies 
{7r(cr)} = limfc-^oo/o-ifeCA) = linife^oo f u \k(A) = {tt(u>)}. Let K be the least integer 
such that uj k+ i ^ a K+1 . Then f^ K = f a]K and f UJK+1 {A) n / (Tif+1 (A) = 0. Since J" 
is injective, each / £ T is injective, which implies that : X — > X is injective. 
Since / w[K = f a \ K , it now follows that f u \ K+ i(A) n = L\K(fu K+1 (A)) H 
L\k(.U k+1 (A)) = L\ K (fa K+1 {A))nf u \ K (f ux+1 (A)) = 0. Since {7r(a)} C / CT |x+i(A) 
and {it(uj)} C /uj|jc+i(j4) it now follows that {7r(er)} n {7r(w)} = 0. 

(v) If SI C T°° is closed then it is compact, because 7 00 is compact. It follows that 
r _1 = 7r|n : i7 — >■ A is a continuous bijective mapping from a compact space SI onto 
a Hausdorff space A. By [T7l Theorem 5.6, p. 167] it follows that t _1 : SI — > A is a 
homeomorphism. It follows that r : A — ¥ SI is a homeomorphism. 

(vi) Suppose that A is connected and A is not a singleton. It follows that SI is 
not a singleton. It follows that SI is not connected. (Since SI contains more than one 
point and is a subset of 7°°, which is totally disconnected when it contains more than 
one point, it follows that SI is not connected.) 

Now suppose t : A — y SI is continuous. Then r is a homeomorphism. It follows 
that A is not connected. But A is connected. So r : A — > SI is not continuous. □ 

In Definition |3 . 3| we define a type of transformation between attractors of iterated 
function systems by composing sections of coding maps with coding maps. We call 
these transformations "fractal" because they can be very rough; specific examples 
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demonstrate that the graphs of these transformations, between compact manifolds, 
can possess a non-integer Hausdorff-Bcsicovitch dimension. 

Definition 3.3. Let J- = (X; /i, / 2 , /at) be a point-fibred iterated function 
system over a compact Hausdorff space X. Let Ajr c X be the attractor of J- . Let 
n jr : I°° — > Ajr be the coding map of T . Let Tjr : Ajr Qjr c 1 00 &e a branch ofnjr. 
Let Q = (Y; gi, <72j <7iv) o point-fibred iterated function system over a compact 
Hausdorff space Y. Let Ag be the attractor of Q . Let ng : L°° — > Ag be the coding map 
of Q . The corresponding fractal transformation is defined to be 

Tjrg : Ajr -> Ag, X I > TTg O Tjr(x). 

In Theorem |3.4| we describe some continuity properties of fractal transformations. 
These properties make fractal transformations interesting for applications to digital 
imaging. 

Theorem 3.4. Let J- and Q be point-fibred iterated function systems as in Defi- 



nition 3.3 Let Tjrg — irgoTjr : Ajr — >• Ag be the corresponding fractal transformation. 

(i) If, whenever o~,u> € fljr, njr(a) — njr(u) ng(a) — ng(uj), then Tjrg : Ajr — > 
Ag is continuous. 

(ii) If Qg '■— Qjr is an address space for Q, and if, whenever er, uj G Qjr, njr(a) = 
ttjt(uj) ^g(cr) — ng(uj), then Tjrg is a homeomorphism and Tj7g = Tgjr. 

Proof, (i) Assume that, whenever a, u) G ftjr, njr(a) = njr(uj) => ng(a) = ng(u>). 
We begin by showing that ngOTjronjr]^ : Cljr — > Ag is the same as irg \-^: : fljr — > Ag. 
Let a G fijr. Then there is u) G fljr such that njr(a) — njr(u:) because Qjr is a 



3.2 



code space for T . Hence Tjr o njr\-rj-~;(cr) = Tjr o 7Tjr|o^-(w) = u> (using Theorem 
(hi)). Hence irg o Tjr o 7rjr|f^r(cr) = ng{uj) = irg (a), where we have used our initial 
assumption. It follows that ng o Tjr o njr\^- = TTg\op- It follows that irg o Tjr o 7Tjr|jy^ 
is a continuous map from a compact space ttjr to a compact Hausdorff space Ag. 
But ng o Tjr o ttjt\-(J— — (ng o Tjr) o jt\q— is the composition of a continuous mapping 

^•HsV • ~~ * ^jr, n ' om a compact Hausdorff space Ojr onto a Hausdorff space Ajr, 
with a mapping ng o Tjr : Ajr — > from Ajr into a Hausdorff space Ag . It follows 
by a well-known theorem in topology, see for example [161 Proposition 7.4, p. 195], 
that Tjrg = irg o Tjr : Ajr — > Ag is continuous. 

(ii) Assume that fig = £ljr is an address space for Q, and that, whenever er, uj G 
fijr, 7rjr(cr) = 7Tjr(ai) <S=> 7rg (er) = 7rg(cj). Then by (i) both of the mappings Tjrg = 
ng o Tjr : Ajr Ag and Tgjr = njrOTg : Ag — > Ajr are continuous. Using the fact that 
the range of Tg is fijr it is readily checked that Tjrg o Tgjr = ng o rjro 7rjr o rg = 
and Tgjr o Tjrg = 7rjr o rgo ng o Tjr — Hence Tjrg is a homeomorphism and 

Zfc - (Tjrg)-\ □ 

Remark 3.5. Suppose that the equivalence relations ~jr and ~gon I°° , induced 
by njr : I°° — > Ag and ng : I°° — > Ag respectively, are the same. Then it is well 
known that, using compactness, the quotient topological space I°° / ^jr— I°° / ^g is 
homeomorphic to both Ajr and Ag . 

See [T] for discussion of relationships between the topology of an attractor and 
the equivalence class structure induced by a coding map. 

4. Construction of sections of n. In order to construct a fractal transforma- 
tion we need to specify a section of n. In order to construct a section of n we have 
to construct an address space f2; that is, we need to specify one element from each of 
the sets in the collection {n^ 1 (x) : x G A}. To do this in a general way seems to be 
difficult; for example, if /i(A) n /2(A) contains a nonempty open set O, then n~ 1 (x) 
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is non-denumerable for all x £ O. However there are two particular related methods. 
These methods yield interesting structure; for example, in both cases the resulting 
address space ft C I°° is mapped into itself by the shift operator, see Theorem |4.5| 
(ii). They are as follows. 

(a) (Masked iterated function system method.) This method requires that J- is 
injective. Define a dynamical system on T : Ajr — > Ajr with the aid of inverses of 
the functions in T . Follow orbits of T to define Tjr : Ajr — > iljr; in effect one uses a 
Markov partition associated with T to define fi.F. 

(b) (Fractal tops method.) Use the dictionary order relation on I°° to select 
a unique element of 7r~ 1 (a;) for each x £ Ajr. This method applies when T is not 
required to be injective. When T is injective, it is a special case of (a). To date, the 
fractal tops method seems to be the easiest to convert to computational algorithms 
and applications. 



4.1. (a) The masked iterated function system method. Definition 4.1 
introduces a special partition of an attractor. 

Definition 4.1. Let J- be a point-fibred iterated function system on a compact 
Hausdorff space. Let A be the attractor of T. A finite sequence of sets M. := {Mi C 
A\i £ 1} is called a mask for J- if 

1. M, C fi(A), i £ I; 

2. M i nM,-=0 ; i,j£l,i^j; 

3. U ieI Mi = A. 

Note that for any x £ A there exists a unique i £ I such that x £ Mi C fi(A). 
This enable us, in Definition |4.2[ to define an associated dynamical system on the 
attractor. 

Definition 4.2. Let {Mi : i £ 7} be a mask for an injective point-fibred iterated 
function system T with attractor A. The associated masked dynamical system for 
J- is 



T:A-> A, { 



(fr\x), x£Mi, 
fz\x), x£M 2 , 



J N \x), x£M 



N- 



Theorem 4.3 associates a unique section of n with a masked dynamical system. 

Theorem 4.3. Let T be an injective point-fibred iterated function system with 
attractor A. Let T : A — > A be a masked dynamical system for J- ' , associated with 
mask M. = {Mi : i £ L}. Let x £ A and let {x„}^L be the orbit of x under T ; that 
is, xq = x and x n — T n (xo) for n = 1, 2, .... Let o~k(x) £ L be the unique symbol such 
that Xk-i £ M ak , for k = 1,2, 3, ... . Then 

= W £ I°°\cr ■= o-(x) = ai(x)a 2 (x)a 3 (x)... £ I°°,x £ A} 

is an address space for J 7 . 

Proof. Let x £ A. We begin by proving that, for all K > 1 and j = 1, 2, K — 1, 
we have 

XK-j £ fa K - J + 1 (x)fv K - J + 1 (x) ° ••• ° f<j K (x){A). (4.1) 

Fix K. We use induction on j. Since Xk-1 £ M aK it follows that xk — f^{ x K-i) so 
xk-i — fa K (x)( x K) £ f<r K {A). It follows that equation 4.1 is true for j = 1. Suppose 
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that equation 4.1 is true for j = 1, 2, J < K— 2. It follows that xk-j £ fcr K - j+i(x)° 
fa K - J+1 (x)°--^7a K (x)(A). We also have xk-j £ M (TR ._ 7+1 so^j = /^jiK-j-i) 

which implies = U k _j(x K -j) £ iW_ JW (f<r K _ J+1 ( x ) ° k-_ J+1 (x) ° ■•• ° /^(x) 

Hence equation |4.1|is true for j = J + 1. This completes the induction on j. 



It follows that 

X = Xq £ f ai (x) ° fa 2 (x)°-° Ia K (x){A) 

for all K. It follows that x £ n^ 1 / CT1 ( a )o/ CT2 ( a ,)o...o/ (7) , ( x )(j4) = 7r(o(x)). It follows that 
x = 7r(o(x)) for all x g A. It follows that (£Im) = that is, tt\q m : r2>i — > A 
is surjective. 

To show that Tt\a M '■ A is injective, suppose o(x) ^ o(y) € f° r some 

x,y £ A. Then for some we have a(x)k ^ which implies M a ^ k n M a ^ k = 0. 

Hence T fc (x) 7^ T fc (y). It follows that x^y.U 

Definition 4.4. Let T be an injective point-fibred iterated function system. The 



address space Qm C 1 00 provided by Theorem 4-3 is called a masked address space 
for J- ' . The corresponding section of 'it , say r : A — ¥ Qm, is called a masked section 
o/tt. 

Masked address spaces and masked sections have all of the properies of address 
spaces and sections, such as those in Theorem |3.2[ and associated fractal transforma- 



tions have the properties in Theorem 3.4 But these objects have additional properties 
that derive from the existence and structure of the masked dynamical system. Some 
of these additional properties are described in Theorem |4.5| 

Theorem 4.5. If T — (X; fx, /2, /jv) is an injective point-fibred iterated 
function system, with attractor A, code space I°° , coding map -k : I°° — > A, mask 
M. = {Mi : i £ I} 7 masked address space VLm an d masked section r : A — > Qm, then 

(i) if T is open then r : A — > is continuous at x £ A if and only ifT k ~ 1 (x) £ 
Int A (M T{x)k ) for all k = 1,2,... ; 

(ii) the shift map 

S : VL M — > f2jK, <7xCT2(73... !-> 02030-4... 

is well-defined, with S(Qm) C 

(Hi) the following diagram commutes 

A ^ A 
r| It (4.2) 

(iv) if there is i £ I such that Mi = fi{A) then S(CIm) = 

Proof, (i) Suppose that I* -1 (a;) £ Int A (M T(x)k ) for all k = 1,2, ... . Let {jW} 
converge to x. Let t(x) = a = 01O2O3... and let t (x^) — = a[ o\ "Vg" .... 
Since x £ Int/SM-o\) there is an integer n\ such that x^ n ' £ M ai for all n > n\. 
Since, for all i £ I, /,~ is continuous (because T is invertible and open) we have 
{T(x(™)) = f„ 1 (x^)}^_ ni converges to T(x) = / CTl (x) as n — ?> 00. Similarly, for any 
given X, there is uk such that T p (x^ n ') £ M a for all p < K and all n > nj(. It 
follows that, for any given i^T, Op = a v for all p < K and all n > nx- It follows 
that {°~^Yn-i conver S es to o, i.e. {r (a^™)) }^ 1 converges to t(x). It follows that 
r : A — > fi^vi is continuous at x. 
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To prove the converse we assume that it is not true that T k ^ 1 {x) € IntA{M T ( x ) k ) 
for all k = 1,2,... . It follows that there is some K > such that T K (x) £ 
JI4 K+1 \M^(M ffK ). Here as in the first part of the proof, we write t(x) = a = 
o~\o~io 3.... It follows that there is a sequence {yj}^ =l that converges to T K (x) with 
Uj M aK+1 for all j. (Any neighborhood of T K (x) must contain a point that is 
in A\M aK+1 .) It follows that {f ai o / ff2 o ... o /^(y.,)}^ converges to f ai o / CT2 o 
••■ fa K (T K (x)) since / CTl o / ff2 o ... o / ffJf : A — > A is continuous. (We define 
fax fa 2 ° ••■ ° /o- = u-) But o / CT2 o ... o (T x (a;)) = x, while t(/ CTi o f aa o ... o 
f<r K {yj))K+i ^ o-K+i because T^/^ o f c% o ... o f aK {yj)) = y 3 £ M„ K+1 . It follows 
that r : A — > fix is not continuous. The desired conclusion follows at once. 

(ii)&(iii) Let O\0~20~-i--- G Let x — t~ 1 (cri<T2<73...). Then roTor -1 ((Tito's...) = 
C2C3C4-", whence <72<73<74... € fix- It follows that the shift map 

S : fix — > Qm7 ~1°~20'3--- !-> CT2Cr3CT4... 

is well-defined, with S(fix) C fix- 

(iv) Let (Tia2<T3... <G fix and let i £ I be such that M, = fi(A). Then A D 
T(A) D T(Mj) = A whence A = T(A). By (in) we have S o r(A) = r o T(A) so 
S(SI m )=t{A) = fix- □ 

We remark that masked dynamical systems are related to Markov partitions in 
the theory of dynamical systems. See for example O Proposition 18.7.8, p. 595]. 

In general a masked dynamical system T : A — >• A depends on the mask Ai. By 
suitable choice of mask we can sometimes obtain a dynamical system with a desired 
feature such as continuity, or which relates the iterated function system to a known 
dynamical system, as illustrated in the following example. 

Example 4.6. Consider the IFS 

T = {M, f (x) = tx, h{x) = -tx + 1} 

where t £ [|, 1) is a parameter. The attractor of A = [0, 1]. Let M\ = [0, 5], M% = 
( 1 , 1] . Then the masked transformation 



T{x) 



x £ Mi, 
x £ M 2 



is continuous. This is the well-known one-parameter tent map dynamical system, see 
for example fTR Exercise 2.4-1, p.78]. We note that, for any x £ (0, 1) the re exist 
a positive integer n such that T k (x) £ [ 0*2 i Jj] f or all k > n (see Fig. As a 



consequence, if A denotes the set of masked addresses of points in [ 2 L 2 1 ; 27 J then the 
masked address space for IF is 

{lVL.Aa\cr £ A,n = 0, 1,2, ...} U {2 111...1ct|c7 eA,n = 0, 1,2,...}. 

n times n times 

Theorem |4. 7| concerns the relationship between masked address spaces correspond- 
ing to distinct masks. It has an application to packing multiple images into a single 
image, as illustrated in Example |5.7| 

Theorem 4.7. Let J 7 be a point-fibred iterated function system with attractor A. 
Let T : A — > A be a masked dynamical system for T corresponding to mask Ai . Let 
[x be a measure on A. Let M! be a mask for J- such that 

M j* G A : {T k (x)}Z f| (ij (M^M'M = j = 0. 
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M 



x 



4> 



Fig. 4.1. Masked dynamical system on the real interva l [0, 1] . The orbit of any point eventually 
enter a trapping interval, [(2t — l)/2t 2 , l/2t]. See Example\4-6\ 




Let tm and tm' be the masked sections of tt corresponding to M and Ai' . Then 
tm{ x ) ^ i~M'(, x ) f or almost all x £ A. 

Proof. If tm{x) — tm'{x) it follows that T k {x) = TJ^cc) for all k, where Tm' '■ 
A — >• A is the masked dynamical system corresponding to M.' . Then T^ A ,(x) ^ 
Uig/ (MjAM-) for all k. But this is impossible for \i— almost all x £ A. □ 

4.2. (b) The fractal tops method. "Fractal tops" is the name we use to refer 
to the mathematics of tops functions, tops code spaces, tops dynamical systems, and 
associated fractal transformations; see for example [SJ [7]- Here we show that, in 
the case where J- is injective, fractal tops arise as a special case of masked iterated 
function systems. Specifically, a tops code space is a special case of masked address 
space, a tops function is a special case of a section of ir, and a tops dynamical system 
is a special case of a masked dynamical system. The computations associated with 
fractal tops tend to be less complicated than those for masked systems. 

Define a dictionary ordering (see for example [T71 p. 26]) on 7°° as follows: if 
ct,ll! £ I°° , a 7^ io, then 



where k is the least index for which ^ ujk, for all er, lo £ I°° . With this ordering, 
every subset of I°° possesses a greatest lower bound and a least upper bound. Since 
tt : I°° — > A C X is continuous and I°° is compact, n^ 1 (x) = {a £ I°° : tt(ct) = x} 
possesses a unique largest element, max7r -1 (izi), for each x £ A. This allows us to 
define an address space C I°° for T by 




a < uj iff cifc > tu k , 



fl = {max7r £ A}. 



The corresponding section of n is 



r : A — > Q, x >— > max7r 1 (x). 



If J- is injective then the tops dynamical system is a masked dynamical system, 
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corresponding to the mask defined by 



M i = f i (A)\\Jf j (A), iel. 
j=i 

The fact that r (x) can be computed from the set tt~ 1 (x) without reference to other 
points on the orbit of x simplifies the computation of r (x) in applications, see for 
example [5,6. Note that we can use the orbits of a tops dynamical system to calculate 
the top address of any point x G A according to t(x) — o\02--- where 

a k =min{ne {1,2,...,JV} : T 3 ^ {x) G f n (A)}. 

5. Applications and examples. 

5.1. Application to image synthesis. Here we generalize the technique of 
color-stealing, introduced in [2] and implemented for example in [181 p. 65-66]. 
Define a picture to be a function of the form 

where £ is a color space. The set D is the domain of the picture. We are concerned 
with situations where D is a subset of an attractor of an iterated function system. 

Let T be an injective iterated function system with attractor A C I 2 . Let T' 
be an iterated function system with attractor A' C ffi 2 and the same code space I°° 
as for T . Let M. be a mask for T and let t_m be the corresponding section of 7r -1 . 
Let 7r' be the coding map for J 7 '. Then we can define a mapping <I> from the space of 
pictures on A' into the space of pictures on A according to 

$mW) =q3WoTM. 

We refer to this procedure as color-stealing because colors from the picture <P' are 
mapped onto the attractor A to define the new picture *P = 



It follows from Theorem 4.5 (i) that the transformation n' o tm is continuous at 
all points x G A whose orbits lie in A\ U^ D _ ri J rk (Jj^id Mj) , where dMi denotes the 
boundary of Mj. In some cases, such as those in Example 5.1 U'^ =Q J rk (Uig/9Mj) is a 



In some cases, such as those in Example 
set of Lebesgue measure zero, so the transformation is continuous almost everywhere. 
This explains patches of similar colors tends to exist in pictures that are obtained 
by color-stealing from real world photos, where patches of similar colors occur for 
physical reasons. 

Example 5.1. Figure 5.1 illustrates color- stealing using (i) fractal tops (left), 
and (ii) a masked iterated function system (right) that is not a tops system. The 
picture 58' from which colors are stolen is Lena embedded in a black surround, shown 
in the middle panel. In (i) T' is an affine iterated function system whose attractor A' 
is a filled square, the domain of 03', such that {fi(A')}f =1 is a set of tiles that tile A' 
by rectangles. In (i) T is the projective iterated function system (KP ; /i, f%, f$, f^), 
where 

,a n x + b n y + c n d n x - 



fn{x,y) = ( 



e n y - 



g n x + h n y + j n g n x + h n y + j„ 



for n = 1,2, 3, 4; 



n 


a n 


K 




d n 




k n 


9n 


K 


jn 


1 


19.05 


0.72 


1.86 


-0.15 


16.9 


-0.28 


5.63 


2.01 


20.0 


2 


0.2 


4.4 


7.5 


-0.3 


-4.4 


-10.4 


0.2 


8.8 


15.4 


3 


96.5 


35.2 


5.8 


-131.4 


-6.5 


19.1 


134.8 


30.7 


7.5 


4 


-32.5 


5.81 


-2.9 


122.9 


-0.1 


-19.9 


-128.1 


-24.3 


-5.8 
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Fig. 5.1. An example of color- stealing using a masked iterated function system (right) and a 
tops function (left). Although the stolen pictures exhibit many discontinuites, the transformations, 
from the original picture of Lena to the stolen pictures, are continuous almost everywhere. See\5. 1\ 



and J- 1 is an affine IFS whose attractor is a filled square; see also J5[ p. 2]. In (ii) the 
iterated function system T and mask M. are the same as in Example \5.5\ while T 1 is 
a perturbed version of Q in Example \5., 



5.2. Fractal homeomorphisms for image beautiflcation, roughening, and 



special effects. Under the conditions of Theorem 3.4 (ii) the fractal transformation 
Tjrg is a homeomorphism. Such homeomorphisms can be applied to pictures to yield 
new pictures that have the same topological properties as the original. For example 
the connectivity properties of the set defined by a particular color will be preserved, 
as will be the property that certain colors lie in an arbitrary neighborhood of a point. 
But geometrical properties, such as Hausdorff dimension and collinearity, may not be 
preserved and, indeed, may be significantly changed. 

Techniques for constructing and computing fractal homeomorphisms using fractal 
tops, with projective, affine, and bilinear IFSs, have been discussed in [SI 151 Rfl 171 IT%] . 
Families of fractal homeomorphisms, built from such transformations in ]R 2 , may be 
established by using code space arguments. Typically, the attractors of the iterated 
function systems in question are non-overlapping, which simplifies the proofs: in some 
situations one only needs to show that the equivalence classes of addresses agree on 
certain straight line segments. The resulting families of transformations are described 
by a finite sets of real parameters. These parameters may be adjusted to achieve 
desired effects such as increased roughness, or continuous (but non-differentiable) 
transformation from a meaningless picture into a meaningful one. Both of these 
effects are illustrated in Example |5.2| 



Example 5.2. Figure 5.3 illustrates three fractal homeomorphisms of the unit 
square applied to Lena. All the iterated function systems involved are constructed 
using bilinear functions defined as follows. Let □ = [0, l] 2 C R 2 denote the unit 
square, with vertices A = (0,0), B = (1,0), C = (1,1), D = (0,1). Let P,Q,R,S 
denote, in cyclic order, the successive vertices of a possibly degenerate quadrilateral. 
We uniquely define a bilinear function B : 1Z — > 1Z such that B(ABCD) = PQRS by 

B(x, y)=P + x(Q-P)+ y(S - P) + xy(R + P-Q-S). 

This transformation acts affinely on any straight line that is parallel to either the x- 
axis or the y-axis. For example, ifB\AB ■ AB — > PQ is the restriction to AB ofB and 
if Q : M 2 -> R 2 is the affine function defined by Q(x, y) = P + x(Q - P) + y(S - P), 
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Fig. 5.2. The four quadrilaterals IEAH, IEBF, IGCD, IGDH, define four bilinear transfor- 
mations and a corresponding iterated function system whose attractor is a filled square. Provided 
that the quadrilaterals are close enough to quadrants of the square, the iterated function systems is 
point-fibred. 




Fig. 5.3. Three fractal homeomorphisms applied to Lena. See Example\5.S\ 



then Q\ab — B\ab- Sufficient conditions for a bilinear iterated function system to be 



point-fibred are given in F^/. Each homeomorphism in Figure 5.3 is generated using 
a pair of iterated function systems of the form in Figure \ 5. <?[ • each such pair has the 
same address structure. 

The following theorem provides practical sufficient conditions for a bilinear IFS 
to be hyperbolic. 

It is more difficult to establish conditions under which pairs of masked overlapping 
iterated function systems, built from geometrical functions such as affines on K 2 , 
yield homeomorphisms. In order to establish that a pair of masked iterated function 
systems provides a fractal transformation that is a homeomorphism, it is necessary 
to establish that the two masked address spaces agree. Interesting non-trivial cases 
involve overlapping iterated function systems, and cannot be reformulated as fractal 
tops. A beautiful family of such examples is provided by Theorem |5.3| 

Theorem 5.3. Let T = ([0,1], f\{x) = ax, f 2 (x) — by + (l — b)), where a > b > 0, 
and a + b > 1, have mask = {M 1} M 2 } where M t = [0,p], M 2 = 0,1], and 

pG [1 — 6, a] . Let Q = ([0,1], giix) = bx, g 2 {x) — ay+(l— a)) have mask M^ 1-23 ) . Then 
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there exists p* E [1 — b, a] such that the corresponding fractal transformation 7>g : 
[0, 1] — > [0, 1] is a homeomorphism when p = p*. The inverse of this homeomorphism 
is Tgjr : [0, 1] — ¥ [0, 1] defined by associating the mask M^~ p *^ with Q. 

Proof. This is an outline. Let D = [0, lfcK 2 . Define affine transformations by 

W- : R 2 -¥ R 2 ,(x,y) ^ (ax, by), 

W + : R 2 -> R 2 ,{x,y) ^ (bx + 1 - b, ay + 1 -y). 

Let S_ = {(x,y) E R 2 : x + y < 1} and S+ = R 2 \S_. We consider the dynamical 
system 



Q : R 2 -> R 2 ,{x,y) ^ 



Wl\x,y) if {x,y)eS-, 
W+\x,y) if (x,y)€S+. 



This possesses a "repeller", a compact set TZ C □, such that 

TZ = Q(TZ). 

In order to define 72., we define 

W : K (□) -> K (□) , C ^ (5_ n W_(C)) U (5^ n W + (C)), 

and 



72= lim ^ fe (D). 

72 is well-defined because it is the intersection of a decreasing sequence of nonempty 
compact sets. (It is quite easy to see that 72 is the graph of a monotone function from 
[0, 1] onto [0, 1].) Using symmetry about the line x + y = 1 and the contractivity of 
W- and W+ in both the x and y directions it can be proved that 72 has the following 
properties. 

(i) Q(72) = 72; 

(ii) 72 is symmetrical about the line x + y = 1 ; 

(iii) P-(72) = P+(72) = [0, 1], where P_ : K 2 -> R denotes the projection in the x 
direction and P + : R 2 — > R denotes projection in the y direction; 

(iv) there is a continuous, monotone strictly increasing function $ : [0, 1] — > [0, 1], 
such that $(0) = 0, $(1) = 1, $(1 - $(x)) = 1 - x for all x € [0,1], and 72 = 
{(i,$(i)):i6[0,l]}; 

(v) there is a unique p* E [1 — b, a] such that <&(p*) = 1 — p*\ 

(vi) there is a continuous, monotone strictly increasing function, ^ : [0, 1] — > [0, 1], 
such that *(0) = 0, *(1) = 1, *(1 - *(y)) = 1 for all y E [0, 1], and 72 = {(y, 1 - 

:ye [0,1]}; 

(vii) *(l-p*) =p*; 

(viii) *(y) = $ _1 (y) for all y e [0, 1]; 

(ix) if p = p*, then the masked dynamical system Tjr : [0, 1] — > [0, 1] obeys 

Tjr{x) = P-(Q(x,<f>(x))) for all x E [0, 1]; 

(x) if p = p*, then the masked dynamical system Tg : [0, 1] — > [0, 1] obeys 



Tg(y) = P+(Q(*(y),y)) for all y E [0,1]. 
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Fig. 5.4. This shows approximations, in shades of grey, to t he " repeller" 1Z of the dynamical 
system Q : R 2 — > M 2 described in the outline proof of Theorem |5.3| An escape-time algorithm, 
similar to the ones discussed in {3) Ch. 7], was used to make this image, with a = 2/3 and b = 1/2. 



These statements imply the theorem. □ 

Remark 5.4. Clearly essentially the same result and proof applies for any anal- 
ogous pair of overlapping strictly increasing functions on [0, 1]. 

Figure [574] illustrates the "repeller". It is a subset of the attractor of the iterated 
function system (□; (ax, by), (bx + 1 — b, ay + 1 — a)), and may be used to compute p* 
as illustrated in Figure [575] 

An example of a fractal transformation, arising from a masked pair of overlapping 



affine iterated function systems, is given in Example 5.5 The resulting homeomor- 
phism with a = | and b = \ yields a picture of Lena with extra large eyes, Figure 

o 

Example 5.5. Let ri p , q :— (X;hi,h2,fi3,h4,) be the family of affine iterated 
function systems defined by X= {(x,y) £ IR 2 |0 < x, y < 1}, r, s £ (0,1), 

h\(x,y) = (rx,ry),h 2 (x,y) = (sx + 1 - s,ry), 

h 3(x, V) = (sx + 1 - s, sy + 1 - s), h 4 (x, y) = (rx, sy + 1 - s). 

Let J- — Ti.2 1 and, for p G [\, |], let M p be the mask{Afi, M2, M3, AI4} where 

Mi = {(x,y) eX :x<p,y<p},M 2 = {(x,y) € X : x > p,y < p}, 
M 3 = {(x,y) el :x>p,y>p},M 4 = {(x,y) e X : x < p,y > p}. 



Let G 



Ui 2. 

2 ' 3 



Then, by Theorem 5.3, there is a value of p = p* € (^, |), such that 

p* then Tjrg : X — > X then is a 



the condition in Theorem 3.4 (U) holds, and if p 
homeomorphism. The value of p* is 

1 2 

p* = max{p e (-, -) : (p, 1 - p) e A*} 

where A* C X is the attractor of TL := (X;(^x,\y),(\x + \, |y 4- |)) illustrated 
in Figure 5.5 Experimentally we find p* = 0.618 which is used to compute Figure 
\1.1\ The original Lena was overlayed on a black background, as in Figure \5.1\ The 
transformed picture also had a black background that has been omitted here. 
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5.3. Application to image filtering. Here we restrict attention to iterated 
function systems defined on X — [0, l] 2 C R 2 . We are concerned with fractal trans- 
formations from [0, l] 2 to itself, associated with a pair of iterated function systems of 
the form 

7- = ([o,i] 2 ,/ 1 ,/ 2) ...,/ w ), g = ([o,i] 2 , 5l , ff2 ,..., ffw ). 

We suppose that Ajr = Ag = [0, l] 2 . In applications to digital imaging, [0, l] 2 is 
discretized. Here we refer to the the locations of pixels as discretized coordinates. We 
suppose that discretized versions of Ayr and Ag have resolutions rjr and rg respec- 
tively. 

Let Pjr : Ajr — > Ajr be a projection operator, namely a function with the property 
Pjr o Pjr(x) = Pjr(x) for all x £ Ajr. For example, 

Pjr(x) = x discretized to resolution rjr, for all x € [0, l] 2 . 

If Tjrg : Ajr — > Ag is a homeomorphism between the non-discretized spaces then 

Pg := Tjrg o Pjr o Tgjr I Ag -> Ag 

is also a projection operator because Tjrg o Tgjr = ig, the identity on Ag, whence 

PG ° Pg = TfG ° P? Tgjr O Tjrg oPjrO Tgjr = Tjrg O Pjr O Ijr O Pjr O Tgjr = ig . 

A trivial but instructive example is provided by choosing T and Q to be the same, 
with TV = 4 and 

h{x,y) = (0.5x,0.5y)J 2 {x,y) = (0.5a + 0.5, 0.5y), 

/sOz, y) = (0.5x + 0.5, 0.5y + 0.5), h(x, y) = (0.5x, 0.5y + 0.5). 

If rjr = rg/2 then Pg is the filter that corresponds to downsampling followed by 
doubling the width and height of each pixel. 

Example 5.6. Let X= {(x,y) £ R 2 \0 < x,y < 1}, p e (0, 1), q = l-p, and 

hi(x,y) = (px,py),h 2 {x,y) = {qx+p,py), 
hz{x,y) = (qx+p,qy + p),h 4 (x,y) = (px,qy+p). 
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Fig. 5.6. Lena before (left) and after (middle) a fractal homeomorphism has been applied. The 
image on the right is the result of applying the corresponding fractal filter. See Example\5.6\ 



The family of IFSs H p := (X; hi, h 2 , /13, hi) has attractor X and address structure 
that is independent of p 6 (0, 1). If we set T = H0.5 and Q = Hq.s then Tj=q : X~^X is 
a homeomorphism. The result of applying the fractal homeomorphism Tjrg to a digital 
(512 x 512,) picture of Lena (left) is illustrated in the middle panel of Figure 5.6 In 
effect the middle image is obtained by composing a projection Pjr t onto a 512 x 512 
pixel grid, with Tjrg . The image on the right is the result of applying Tgjr to the 
middle image, that is, the right-hand image is the result of applying the projection 
operator, that we may refer to as a "fractal filter" , Tgjr o Pjro Tjrg to the original 
Lena. 

5.4. Application to packing multiple images into a single image. Sup- 
pose we have a collection of masks {Mj} on A such that the conditions of Lemma 
4.7 are holds true for any pair of masks. Let the second iterated function system T' 
be such that almost all points of the attractor A' have a unique address. Thus by 
Lemma 14.71 



vr' o t m .(x) ^ tt' o T Mk (x) 

for almost all x € A and j ^ k. Therefore we can consider the collection of pictures 
&Mj(%$') as an almost disjoint fractal decomposition of the picture *B'. 

This leads us to the following trick. By means of fractal transformations we map 
different pictures to different almost disjoint components of 05'. By inverting the 
transformations we retrieve (approximations to) the original pictures, as illustrated 
in the following example. 

Example 5.7. See Figure \5~1\ Let 



T = (□, fi(x, y) = (y, 0.6(1 - x)), f 2 (x, y) = (y, 0.4 + 0.6a;)) 



and 



g = (□, h(x, y) = (y, 0.5(1 - x)), f 2 (x, y) = (y, 0.5 + 0.5a:)). 
We define a family of masks for J- by 

M p = {A'h = {(x,y) G □ : x<p},M 2 - {(x,y) £ □ : x>p}}, P e [0.4,0.6]. 



□ is the attractor of both systems. We denote the corresponding masked fractal trans- 
formation Tjrg by Tjfg for p G [0.4,0.6]. Then T^g is injective 
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Fig. 5.7. The pictures on the right and left were obtained by masking; the "encoding" is shown 
in the middle image. In this example each of the image is of resolution approximately 512 X 512. 
See Example\5.' 



and invertible on its range, Tj^g(\3). Moreover, from Theorem J^.l it follows that 
A(Tj?J(D) n T^(D)) = for all p ^ q, where A is Lebesgue measure. Let *}3 and 



be two pictures, each supported on □. Then TfejpP) = CP o (Tj?g) 1 is a picture 

supported on T^(D) and T^(0) = o (T^)" 1 is a picture supported on Tj?g(\3). 
We choose p = 0.44, q — 0.56, ^5 =Lena and =Inverted-Lena, where both pictures 
are 512 x 512. The digitized combined picture := T ^l (ffi ) U T^j(0), also of reso- 
lution 512 x 512 7 is shown in the middle panel of Figure \5. 7| Pixels which correspond 



to points in T^(D) both Ti?i(n) are colored white. The left-hand panel in Figure 



,(</) 



5. 7 illustrates the picture 9{oT^} and the right-hand panel shows tRoT^l . Hence we 



can ''store 7 ' the two pictureswe have that for each distinct choice of p £ [0.4,0.6] the 
ranges of the M0.44 a nd .Mo. 56 respectively. 

5.5. Measure theory image packing. Here we describe a different method 
for encoding several images within a single image. The method strives to create a 
single image that is simultaneously " homeomorphic" , under different fractal transfor- 
mations, to several different given images. An example of such a single "encoded" 
image is shown in the middle panel of Figure |5.8| two different fractal homeomor- 
phisms applied to it (but in the digital realm) yield the images on the right and 
the left. The underlying idea is that different invariant measures, belonging to the 
same iterated function system, but associated with different probabilities, are mutu- 
ally "singular continuous", and are concentrated on different sets, albeit they have 
the same support. The method uses the chaos game algorithm, both to encode and 
decode. 

Consider for example the IFS T = ([0, 1]; h{x) = 0.5a:, f 2 (x) = 0.5x + 0.5). If 
we associate probabilities p\ = 0.5 and p2 — 0.5 with f\ and fi respectively, then 
the associated Markov operator has invariant probability measure fj, equal to uniform 
Lebesgue measure supported on [0, 1], concentrated on the set of points whose binary 
expansions contain equal proportions of zeros and ones. On the other hand, if we 
associate probabilities p\ =0.1 and P2 = 0.9, then the associated Markov operator 
has invariant probability measure fx that is also supported on [0,1], but singular 
continuous and concentrated on the set of points whose binary expansions contain 
nine times as many ones as zeros. Hence, if the chaos game is applied to J- using 
the first set of probabilities, the points of the resulting random orbit will tend to be 
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disjoint from those obtained by of applying the chaos game using the second set of 
probabilities. The precise manner in which the orbits of the two systems concentrate 
is governed by convergence rates associated with the central limit theorem. The same 
idea can be applied to "store" several images in a single image E: a stored image is 
retrieved by applying the appropriate fractal transformation to E. 

We describe the method by means of an example. 

Example 5.8. We use three iterated function systems: 

F = {^,h,h,h,fi), Q = (□,51,52,53,34), U = (0,h 1 ,h 2 ,h 3 ,h 4: ) where 
h{x,y) = (0.66z,0.34y), f 2 (x,y) = (0.34x + 0.66, 0.3%), 
f 3 (x, y) = (0.34a; + 0.66, 0.66y + 0.34), f 4 (x, y) = (0.66a;, 0.66?/ + 0.34); 
gi(x,y) = (0.34a;, 0.66jy), g 2 (x, y) = (0.66a; + 0.34, 0.6%), 
g 3 (x, y) = (0.66a; + 0.34, 0.34y + 0.66), g 4 (x, y) = (0.34a;, 0.34y + 0.66); 
hi(x,y) = (0.5a;, 0.5y), h 2 (x,y) = (0.5a; + 0.5, 0.5y), 
hz{x, y) = (0.5a; + 0.5, 0.5y + 0.5), h 4 (x, y) = (0.5x, 0.5y + 0.5). 

Then T is associated with probabilies P = {pi = area (/i(D)) : i = 1,2, 3, 4} Similarly 
the iterated function system Q is associated with probabilities P = {pi = area(gi(\3)) : 
i = 1,2,3,4}. 

The goal is to "store" two standard digital color images, Pepper and Lena, each 
512 x 512, in a single color image E, also 512 x 512 . The image E is supported on 
□ and associated with two probability measures, and Jl-u, invariant under T-L with 
probabilities P and P respectively. 

(1) In order to "encode " Pepper, we run a coupled chaos game algorithm with 
T associated with Pepper, supported on a copy of □, and T-L associated with E, with 
(supposedly) i.i.d. probabilities P; that is, we compute a random sequence of points 
\{X k ,Z k ) G □ X □ : k = 0,1,2, ...if} where K = 10 6 , where X = (0,0) G □ 
(associated with Pepper) and Z$ = (0,0) € □ (associated with E) and 

X k = f ak (X k ^),Z k = h ak {Z k _ x ) for k = 1,2,..., K, 

where a k —i with probability p t . At each step the pixel containing Z k in the (initially 
blank) image E is plotted in the color of Pepper at the pixel containing X k . At the 
end of this process E consists of an "encoded" version of Pepper. 
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(2) In order to "encode" Lena, we again run a coupled chaos game algorithm with 
Q associated with Lena, supported on a copy o/D, and % associated with E (already 
"painted" with an encoding of Pepper) with probabilities P; that is, we compute a 
sequence of points {(Y k ,Z k ) G □ x □ : k = 0, 1, 2, ...K} where K = 500,000, where 
Yq = (0,0) <E □ (associated with Lena) and Zq = (0,0) G □ (associated with E) and 

Y k = f ak (Yfc-O, Z k = K h (Z k _i) for k = 1, 2, K, 

where a k = i with probability pi. At each step the pixel containing Z k in the image E 
is plotted in the color of Lena at the pixel containing Y k . ( The pixel is overwritten by 
the latest colour value.) We have used half as many iterations in the encoding of Lena 
as we did for Pepper, because a proportion of the points that correspond to Pepper are 
overwritten by points corresponding to Lena. 

An image E that is a realization of steps (1) and (2) is shown in the central panel 



of Figure 5.8 The image E is approximately homeomorphic to both of the images 
Pepper and Lena, under the fractal transformations Tjr n and Tg-u respectively, that 
is 

Pepper = Tj^-h(E) and Lena = Tgu(E). 

In practice, to obtain the decoded images, shown on the left and right hand sides of 
Figure [5TiS[ we use the chaos game algorithm again, as follows. 

(3) In order to "decode" Pepper, we run a coupled chaos game algorithm, with 
probabilities P, with T associated with an image (initially blank), supported on a copy 
of O, and % associated with E (now encoding both Pepper and Lena). That is, we 
compute a random sequence of points {{X k ,Z k ) € □ X □ : k = 0,1,2, ...K} where 
K = 10 6 , where Xq = (0,0) G □ (associated with Pepper) and Zq = (0,0) G □ 
(associated with E) and 

X k = U k (X k -i),Z k = h ah {Z k _ x ) fork = 1,2,.. .,K, 

where a k — i with probability pi. At each step the pixel containing X k in the (initially 
blank) copy o/D is plotted in the color of E at the pixel containing Z k . The result of 
such a decoding, starting from the encoded E illustrated in the middle panel, is shown 
in the left panel in Figure \5.8\ 

(4) In order to decode Lena, we run a coupled chaos game algorithm, with prob- 
abilities P, with Q associated with an image (initially blank but to become the de- 
coded image), supported on a copy of □, and W associated with E. That is, we 
compute a random orbit {{Y k , Z k ) G □ X □ : k = 0, 1, 2, ...K} where K — 10 6 , where 
Yq = (0, 0) G □ (associated with Pepper) and Zq = (0, 0) G □ (associated with E) and 

Y k = f* k (Xk-i),Z k = K k {Z k _ x ) fork = 1,2,..., K, 

where o~ k = i with probability pi. At each step the pixel containing Y k in the (initially 
blank) copy o/D is plotted in the color of E at the pixel containing Z k . The result of 
following this decoding algorithm, starting from the encoded E in the middle panel, is 
shown in the left panel in Figure \5.< 
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